Measuring the characteristic function of the work distribution 
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We propose an interferometric setting for the ancilla-assisted measurement of the characteristic function of 
the work distribution following a time-dependent process experienced by a quantum system. We identify how 
the configuration of the effective interferometer is linked to the symmetries enjoyed by the Hamiltonian ruling 
the process and provide the explicit form of the operations to implement in order to accomplish our task. We 
finally discuss two physical settings, based on hybrid opto-/electro-mechanical devices, where the theoretical 
proposals discussed in our work could find an experimental demonstration. 

PACS numbers: 



Thermodynamics is one of the pillars of physical, chemical 
and biological sciences. Its principles are able to predict and 
explain the occurrence and efficiency of complex chemical re- 
actions and biological processes. In physics, the conduction of 
heat across a medium or the very concept of the arrow of time 
are formulated thermodynamically. In information theory, the 
definitions of information and entropy are also given in ther- 
modynamical terms and explicitly related to this framework. 
Even more, the tightness of the link between information and 
thermodynamics can be deduced from the thermodynamical 
interpretation of the landmark embodied by Landauer's era- 
sure principle 111 • 

The outstanding levels of dexterity and reliability that char- 
acterise the current experimental ability to control and guide 
processes at the microscopic scale opens up a series of tan- 
talising questions, the most pressing one being probably the 
following: what happens to the principles of thermodynam- 
ics when we deal with the quick, non-quasistatic dynamics of 
small quantum systems brought out of equilibrium? An in- 
valuable tool for the formulation of an answer in this sense has 
been provided with the formulation of non-equilibrium fluctu- 
ation relations and their subsequent extension to the quantum 
formalism |2, 3], which has very recently enabled interest- 
ing investigations at the border of quantum physics, thermo- 
dynamics, and information theory [4], including the experi- 
mental proposal of building blocks for quantum thermal ma- 
chines ||5[], the study on the connection between quantum fluc- 
tuation relations and critical phenomena in many-body quan- 
tum systems {6,2}, the verification of key features of the quan- 
tum thermodynamical framework as embodied by the Jarzyn- 
ski equaUty OS MlOll . and the extension to open-system dynam- 
ics Lnl. 

It is precisely the Jarzynski identity to serve the motivations 
behind our work: the key point for its verification and use, 
both in the classical and quantum formulations 11^,^7], is the 
determination of the work distribution following a process un- 
dergone by a system. It is thus central to the potential applica- 
tion and development of the investigation on non-equilibrium 
quantum dynamics. To the best of our knowledge, the most 
advanced experimental scheme for the determination of work 
distribution following a quantum process relies on hard pro- 



jections onto energy eigenstates of the system directly, which 
can be of significant practical difficulty |8]. It is thus impor- 
tant to devise strategies for reaching this outstanding goal. 

Our work moves along these lines to propose a way to in- 
fer the quantum statistics of a work distribution by relying 
on a well-known interferometric approach that delegates the 
retrieval of the information we are after to routine measure- 
ments performed on a finite-size ancilla. In fact, we demon- 
strate how a qubit-assisted Ramsey-like scheme is effective in 
fully determine the characteristic function of the work distri- 
bution following a general quantum process. This is enough 
information to access the work distribution itself and, thus, 
employ the framework of fluctuation relations to an out-of- 
equilibrium configuration. We clearly identify the relation 
between symmetries that the quantum process might enjoy 
and the corresponding Ramsey interferometer Although our 
scheme does not rely on the details of a specific experimental 
setting and can thus be applied to a vast range of interesting 
physical situations (cf. Ref. fl2] for a related analysis applied 
to an ion-trap setting), as an illustration we apply it to the case 
of a (micro-/nano-)mechanical oscillator coupled to an effec- 
tive two-level system and undergoing a process that displaces 
its state in the phase space, thus showing the effectiveness of 
our scheme to a situation of current strong experimental inter- 
est. We believe that designing experimentally viable ways to 
access quantum statistics of non-equilibrium dynamical pro- 
cesses can embody a significant step towards the experimental 
grounding and the development of this fascinating area of in- 
vestigation and the spurring of potential ramifications in areas 
such as quantum control and foundations of quantum mechan- 
ics 

Quantum fluctuation relations: a brief review.- In order to 
provide a self-contained discussion, here we give a brief sum- 
mary of the quantum fluctuation formalism that will be used in 
the rest of the manuscript. We consider a process described by 
a Hamiltonian 'HiAit)) that depends on a so-called work pa- 
rameter /1(f), which is assumed to be externally controlled. A 
system S is assumed in contact with a reservoir and initialised 
in a thermal equilibrium state p'^{A.o) - e^^'^'-'^''^ Z.{Ao) at in- 
verse temperature /3 and work parameter A{Q^) - Ao [here 
Z(A) - Tre"'^'^*''* is the partition function]. At t - 0+, we 
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detach S from the reservoir and perform the process, which 
consists in changing the work parameter to its final value 
A{t) = Aj. It is now convenient to decompose the two Hamil- 
tonians connected by the process as "//(/lo) - 2„ E„{Aq) \n) {n\ 
^nd<ik{Ar) = Y.,nKi^r)\m){m\, where {E„,\n)) [(£;, |m»] 
is the [to*] eigenvalue-eigenstate pair of the initial [fi- 
nal] Hamiltonian. The corresponding work distribution can 
be written as L3J 



(1) 



Here, we have introduced the probability /7° that the system is 
found in state \n) at time f = and the conditional probability 
p^^l^j to find it in state |to) at time f = r if it was in \n) at f = 
and has evolved under the action of the propagator U{t,Q). 
Eq. ([TJ encompasses both the thermal statistics of the initial 
thermal state (given by p^l) and the fluctuations arising from 
quantum measurement statistics (given by p]^^^)- The role of 
P^{W) becomes immediately evident when considerin g th e 
Tasaki-Crooks relation P^(W)/P^(-W) = e^"'-*'') iHIH, 
which states that the ratio between the forward work distri- 
bution P^{W) and the backward one P^(-W), obtained im- 
plementing the process Ar — > Aq, starting from p'^iAj) and 
evolving through W(t,Q), is related to the difference in the 
equilibrium free-energy h.F of the system. For our purposes, it 
is convenient to define the characteristic function of the work 
distribution [16(1 



Tr U\t, 0)e'"^('''>f/(T, 0)e-"'"^'"'VG(^o) 



-iu'HiA{)) 



(2) 



From Eq. the Jarzynski equality JlCl is found by introduc- 
ing the parameter u - ifi 'm such a way that 



-/iAF 



(3) 



Eq. (O demonstrates the central role played by W (and, in 
turn, Pf{W) and its characteristic function) in the determi- 
nation of the equilibrium properties of a system driven by a 
process. In what follows, we shall illustrate a protocol for 
the interferometric determination of xiu, t), which would then 
enable the convenient evaluation of the figures of merit dis- 
cussed above. 

A simple illustrative case.- In order to fix the ideas before at- 
tacking the general protocol and illustrate the idea behind our 
proposal in an easy-to-grasp situation, we consider a Hamil- 
tonian for S having the form 



9{s(t) ^ g(A,)h, 



(4) 



with h being the operatorial part of the Hamiltonian that re- 
mains unchanged irrespective on the process implemented for 
the change of the work parameter (notice that, for easiness of 
notation, from now on we will replace A(t) — > A,) and speci- 
fied fully by the assignment of the function g(A,). As a conse- 
quence, the system Hamiltonian commutes with itself and the 



time-evolution operator U{t,Q) = g-'T'^sW' ^ ^-ih£g(A,)dt 
all instant of time. That is 

[nki,nkf^ = and [U{T,Q),i<i] = [t/(T, 0),'///] = (5) 

with iki = iksiO) = g{Ao)h YHf = iksir) = g{A,)h]. Due to 
such symmetries, it is straightforward to check that the char- 
acteristic function of the work distribution simplifies to 



;rsimp(M) = Tr 



= Tr 



f/"'"(r, 0)e'"''^^(-*^>f/(T, 0)e-'''''^>^^'^p{AQ) 



(6) 



The interpretation of this result is straightforward: the char- 
acteristic function is fully determined by the changes induced 
by the process responsible for the variation of the work pa- 
rameter in the Hamiltonian of the system S. In turn, this 
allows us to make a significant step forward in the illustra- 
tion of the proposed interferometric setting. Indeed, let us 
now introduce a simple qubit system playing the role of an 
ancilla A whose role is to assist in the measurement of;i^sump- 
Moreover, we consider the system-ancilla evolution G{u)V{u), 
where V(u) - e^''^'" (8> 1a is a local transformation on the S 
Hilbert space and G(u) is the controlled AS gate 



G(u) ^Is® \0}a{0\ + e 



-i(rHf-n)ii I 



|1>A<1| 



(7) 



which applies the transformation to the state of 

S only when the ancilla is in |1)a and leaves it unaffected 
otherwise. Inspired by Ramsey-like schemes for parameter 
estimation 11171 llSll . our protocol proceeds as follows: We 
prepare the ancilla in |0)a and apply a Hadamard transform 
H - {&x + a-.)l V2 [24] that changes it into the eigenstate of 
the ;c-PauH matrix |-i->a = (|0)a + |1)a)/ V2- We then apply the 
sequence G{u)V{u) on the system-ancilla state |-i-)a ® p'^ 
and subject the ancilla to a second Hadamard transform [cf. 
Fig. [Tl. The controlled system-ancilla gate G(u) in general 
establishes quantum correlations between the qubit A and S . 
The fingerprint of this is provided by the fact that information 
on the system can be retrieved from the ancilla. In fact, by 
tracing out the state of 5 , we find 

PA = Txs[HG{umu){\+)A <+l ®pt)V\u)GHu)H] 

1 . (8) 

with a - Re;t'simp and v = Im;t'simp. This proves our claim and 
the effectiveness of our protocol for the reconstruction of the 
characteristic function of the work distribution: the tomogra- 
phy of populations and coherences of the ancilla state enables 
the full determination of ;^'simp(M). 

General protocol- We now relax the assumption made previ- 
ously on the form of the Hamiltonian of the system and con- 
sider the general case embodied by the conditions 



VHi, ikf] + and [t/(T, 0), -H.^/,] + 0. 



(9) 
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FIG. 1: (Color online) Quantum circuit illustrating the simplified 
protocol for the measurement of the characteristic function of the 
work distribution through a single-probe interferometer The ancilla 
A is a qubit initialised in 10)^ and undergoes a Hadamard gate H. 
The system S is prepared in a thermal state and is subjected to 
the local transformation V . The controlled gate G, followed by a 
second Hadamard gate on A, implement our protocol. We refer to 
the body of the manuscript for the explicit form of the gates (whose 
dependence on u has been omitted for easiness of notation). 



FIG. 2: (Color online) Quantum circuit illustrating the protocol for 
the measurement of the characteristic function of the work distribu- 
tion through a single-probe interferometer in the most general case. 
In order to ease the visualization of the scheme, here we have ab- 
sorbed the Hadamard gate H that should be used as the start and end 
step of the protocol into the preparation of the ancilla state and the 
basis onto which the measurements should be performed. The sys- 
tem 5 is assumed to be prepared in a thermal state p'!! . The sequence 
GiVGi implements our proposal. 



Correspondingly, the characteristic function takes the full 
form given in Eq. (|2]l. In this case, the general interferomet- 
ric approach illustrated above still applies, the only difference 
being the form of the controlled operation to be applied on the 
state of the system. More explicitly, we now need to imple- 
ment the joint AS gate 

G{u, T) = f/(T, 0)e-'^'" |0>A<0| + e-''^'"lJ{T, 0) |1>^<1| , (10) 

which can be broken down into the concatenation of a local 
transformation on S and two A-controlled gates as G{u, r) - 
G2(u, t)V{u)G\{u, t) with 

Gi(m,t) = is ® |0)^<0| + e'^'"f/(T,0) ® 1 1)^(1 1 , 

y(M) = iA®e"''^'V^''', (11) 

G2{u,t) = tj{T,Q)e>*'" ® 10)^(01 + is® |1>A<1| • 

The sequence of operations to apply to the S -A system is illus- 
trated in Fig.|2] which provides the quantum logic circuit for 
the interferometric protocol addressing the most general pro- 
cess undergone by the system. Using the same preparation of 
the ancilla as above and the two Hadamard transforms, we ob- 
tain a reduced state of A identical to Eq. (|8]l with the replace- 
ments a — > Re;\f(M, t) and v — » l-mx{u, t) with x(m, t) as in 
Eq. (|2]l. An appropriate sequence of measurements performed 
on the ancilla state only is thus able to provide information on 
the work distribution following the process undergone by S . 
Potential physical examples. -Among others [12], two explicit 
situations of enormous current experimental interest can well 
be used in order to illustrate the main findings of our work. 
They are both based on hybrid coupling configurations be- 
tween an effective two-level system and a mechanical oscil- 
lator, which can be either microscopic (in a hybrid cavity op- 
tomechanics setup) or nanoscopic (for an electromechanical 
implementation). In what follows, we briefly show how to 
achieve the effective Hamiltonians regulating the processes 
that we have in mind in both such scenarios and illustrate the 
principles of our proposal by explicitly calculating the char- 
acteristic function following it. 

We start from a microscopic setting, where we consider a 
three-level atom in a A configuration coupled to a single-mode 



optical cavity that is pumped by a laser field at frequency coj, 
and having a movable mirror. The atom is driven by a sec- 
ond field at frequency w, that enters the cavity radially [cf. 
Fig. [3](a)]. The logical states {|0),|1)) of the ancilla system 
A are encoded in the fundamental atomic doublet (|e) being 
the common excited state). The scheme includes the diiving, 
at rate Q, of the transition |l)<-^|e) by the external field at fre- 
quency oji. The transition |0)<->|e) is guided by the cavity field 
(frequency Wc), at rate g. Both the fields are detuned by 6 
from the excited atomic, while the detuning between cavity 
field and external pump is A=Wf-Wp. 

System S is embodied by the movable mirror, which os- 
cillates harmonically at frequency ojs and is driven by the 
cavity field through radiation-pressure (at a rate x) We 
assume large single-photon Raman detuning and negligible 
decay rate from the atomic excited state, so that 6 » 
(O, g) » -ye and an off-resonant two-photon Raman transi- 
tion is realized. We assume A» g,x, so that both the atomic 
excited state and the cavity field are only virtually populated 
and they can be eliminated from the dynamics of the system, 
and move to an interaction picture defined by the operator 
(jjpC^c + ojj \e){e\ + wio |0)yi(0| with (c, c^) the bosonic oper- 
ators of the cavity field. 

We thus obtain the effective Hamiltonian ll20ll (we take 
units such that ft-l throughout the paper) '//micro = cjsb^b + 
A\\)A{\\®{P+b)WVi\\A= Xg'^Qrl5^A^ and (^, S+) the bosonic 
operators of the mechanical oscillator. Through the two- 
photon Raman transition, the virtual quanta resulting from the 
atom-cavity field interaction are transferred (by the bus em- 
bodied by the cavity field) to the mechanical oscillator. The 
state of the latter is coiTespondingly displaced in its phase 
space, in a way that is controlled by the state of A. By driv- 
ing the cavity with a bichromatic pump with frequencies Upit) 
and u)p{t) + a>m and relative phase 0(f), the effective coupling 
between the ancilla and the system becomes such that dis- 
placements in any direction of the phase space of the mov- 
able mirror can be arranged ll2ll - l23l . Moreover, consider- 
ing a time-dependent amplitude of the driving field, we get 
A ^ A, - xg^^^{t)/S^A^, so that the model 

Kicro = '^sPb + ^, ID^dl ® + Se-''^«) (12) 
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is finally obtained. As our examples aim to only have an il- 
lustrative character, we neglect the sources of damping and 
decoherence in the system, whose inclusion will require the 
extension of Eq. ^ to open-system dynamics. While this is 
possible due to the advances reported in 111 ill , it goes signif- 
icant beyond the goals of this paper and will be left for a 
future investigation. However, it is important to stress that 
current progresses in the fabrication of mechanical oscillators 
allow for very small decoherence rates, while optical cavities 
with large quality factors are currently employed in optome- 
chanical experiments |19], thus making a quasi-unitary pic- 
ture plausible indeed. The state of the atomic system can be fi- 
nally manipulated and reconstructed through an optical probe 
and employing standard tools in quantum optics. 

A similar effective model is obtained by considering the 
second physical setting that we have in mind and shown in 
Fig- [3(b). This involves a nanomechanical oscillator (a nana 
beam) coupled capacitively to a Cooper-pair box (CPB) op- 
erating at the so-called charge degeneracy point ll26ll . where 
the dynamics of the CPB can be approximated to that of a 
two-level system encoded in the space spanned by states \a±} 
that are symmetric and antisymmetric superpositions of states 
with exactly and 1 excess Cooper pairs in the large super- 
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FIG. 3: (Color online) (a) Hybrid micro-optomechanical setting for 
the achievement of an effective Hamiltonian suitable for the verifica- 
tion of the protocol for the measurement of the characteristic func- 
tion of the work distribution. The process is undergone by a system S 
embodied by the movable mirror of the optomechanical cavity. The 
ancilla A of our protocol is encoded in the ground-state doublet of a 
three-level atom, (b) Nano-mechanical version of the setup. We con- 
sider 5 as embodied by an electrically driven nanobeam (bias voltage 
Vv)- The ancilla is a CPB coupled to S via the capacitance C, . The 
state of the CPB is controlled by the gate voltage Vg (coupled to the 
box through the capacitance Cg) and the Josephson energy Ej. We 
work at the charge degeneracy point. 



conducting island shown in Fig. [3] (b) and encode our an- 
cilla. The natural Hamiltonian of the system reads "Hi = 
(6 ft(0) (p+(x)sb'^b with Q and (p the canonical charge 

and phase operator of the CPB, C, the total capacitance of the 
island, gg(f) = CgVg{t) + C^^VAt) the total gate charge, Ej the 
Josephson energy and Q the frequency of the nanomechani- 
cal oscillator 12611 . By defining £t — \a+){a-\ + ex- 
panding "Hi in series of the ratio between the actual position 
of the mechanical oscillator and its equilibrium distance from 
the CPB (the amplitude of the oscillations is assumed to be 
small enough that only first-order terms are retained in such 
expansion) and adjusting the gate and driving voltages such 
that Qg{t) ^ 0, the interaction Hamiltonian of the system can 
be cast into the form 



(13) 



(the explicit form of At is inessential for our tasks) |25]. The 
state of the ancilla can be processed by tuning the gate volt- 
age Vg{t) and the Josephson energy, while single-qubit mea- 
surements in this scenario are performed using single-electron 
transistors |26]. 

Both the models derived briefly above entail a conditional 
displacement of the harmonic oscillator Actually, the phase- 
space displacement is precisely the time-dependent process 
that we choose for S . We thus take the coupling strength A, 
from Ao = to Aj at the end of the process. In what fol- 
lows, we concentrate on the model embodied by Eq. ( fTSl l. 
which can be easily diagonalised by the means of appropri- 
ate, time-dependent Weyl displacement operators. A lengthy 
yet straightforward calculation based on standard phase-space 
methods allows us to fully evaluate the density matrix of the 
ancilla associated with the process experienced by the har- 
monic oscillator and following the protocol that we have sug- 
gested, which can be conveniently cast into the form 



;r(M,T) = e-2'«^^^™^^lTrs[e 



iuiLOsb^ b+Ar(b^ +h)) th 



^ (n -H l)«+i " ^ 

(14) 

where X!f'(jc) are associated Laguerre polynomials of order 
(n, p) and argument x and we have introduced the parameters 
yiAu) = ^^(1 +e''^'"')/w5, Jr = A,e-'"^''dt, n = lUe^'^' - 1). 
The details of this derivation, and an explicitly analysis of the 
characteristic function will be given in Ref. ||27ll. It is worth 
noticing that the protocol for the measurement of the charac- 
teristic function performed using the model in Eq. (fTjt is en- 
tirely equivalent to the above, modulo local unitary operations 
performed on the ancilla qubit [27]. 

Conclusions. -'^t have proposed an ancilla-assisted interfer- 
ometric protocol for the measurement of the characteristic 
function of the work distribution corresponding to a process 
that is enforced on a system. The scheme requires both lo- 
cal operations performed on the oscillator and joint gates con- 
trolled by the state of the ancilla, and shares similarities with 
Ramsey-based strategies for parameter estimation. We have 
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illustrated the working principles of our protocol by tack- 
ling the case of a mechanical oscillator undergoing a time- 
dependent displacement in phase space and coupled to an 
atomic ancilla, which embodies an interesting case for non- 
equilibrium quantum dynamics of current strong experimental 
interest. It should be noted, though, that our proposal bears no 
dependence on the specific experimental setting being consid- 
ered and is widely applicable to any physical system allowing 
for a controllable and engineerable system-ancilla interaction 
and easily performable measurements on the latter (cf. lil2ll 
for an example based on an ion-trap setup). 

As the characteristic function of the work distribution is a 
key element in the framework of quantum fluctuation relations 
for out-of-equilibrium system, designing experimentally vi- 
able strategies for its inference is an important step forwards 
in the grounding and development of out-of-equilibrium quan- 
tum thermodynamics. We believe that our proposal con- 
tributes to such a quest in a relevant way and opens up the pos- 
sibility for an experimentally implementable verification of 
the connections between out-of-equilibrium quantum statis- 
tics and criticality in a quantum many -body system ll6i [l3l[l8l l. 
Very interesting routes for the application of our protocol in- 
clude the study of quantum statistical properties of quantum 
thermal machines 11411 . 
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